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ABSTRACT 



Context. The star formation process in large clusters/associations can be strongly influenced by the feedback from high mass stars. 
Whether the resulting net eff'ect of the feedback is predominantly negative (cloud dispersal) or positive (triggering of star formation 
due to cloud compression) is still an open question. 

Aims. The Carina Nebula complex (CNC) represents one of the most massive star-forming regions in our Galaxy. We use our 
Herschel far-infrared observations to study the properties of the clouds over the entire area of the CNC (with a diameter of ^ 3.2°, 
which corresponds to ^ 125 pc at the distance of 2.3 kpc). The good angular resolution (10'' - 36") of the Herschel maps corresponds 
to physical scales of 0.1 - 0.4 pc, and allows us to analyze the small-scale (i.e. clump-size) structures of the clouds. 
Methods. The full extent of the CNC was mapped with PACS and SPIRE in the 70, 160, 250, 350, and 500 yum bands. We determine 
temperatures and column densities at each point in this maps by modeling the observed far-infrared spectral energy distributions. We 
also derive a map showing the strength of the UV radiation field. We investigate the relation between the cloud properties and the 
spatial distribution of the high-mass stars, and compute total cloud masses for diff'erent density thresholds. 

Results. Our Herschel maps resolve, for the first time, the small-scale structure of the dense clouds over the entire spatial extent of 
the CNC. Several particularly interesting regions, including the prominent pillars south of rj Car, are analyzed in detail. We compare 
the cloud masses derived from the Herschel data to previous mass estimates based on sub-mm and molecular line data. Our maps 
also reveal a peculiar "wave"-like pattern in the northern part of the Carina Nebula. Finally, we characterize two prominent cloud 
complexes at the periphery of our Herschel maps, which are probably molecular clouds in the Galactic background. 
Conclusions. We find that the density and temperature structure of the clouds in most parts of the CNC is dominated by the strong 
feedback from the numerous massive stars, rather than random turbulence. Comparing the cloud mass and the star formation rate 
derived for the CNC to other Galactic star forming regions suggests that the CNC is forming stars in an particularly efficient way. We 
suggest this to be a consequence of triggered star formation by radiative cloud compression. 

Key words. ISM: clouds - ISM: structure - Stars: formation - ISM: individual objects: NGC 3372, Gum 31 , G286.4-1 .3, G289. 0-0.3 
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l/^ 1. Introduction millimeter survey of the CNC with LABOCA at the APEX 

telescope revealed the structure of the cold a nd dense clouds in 

CnThe Carina Nebula complex (CNC hereafter) (e.g. the complex in detail (iPreibischetaLllffii). 
^— r ISmith & Brooksl 120081 for an overview of the region) is , ^ ^ • . • . . , . , , , ,- , 

^ bne of the richest and largest high-mass star forming regions in , The CNC is an interesting region in which to study the feed- 

our Galaxy. At a moderate distance of 2.3 kpc, it contains at ^^""^ ^^^""^^ niimerous very massive and luminous stars 

> least 65 O-tvoe stars (Smith 2006) and four Wolf-Rayet stars *e surrounding clouds. The strong lonmng radiation and 
• ^ ^ggg ^ ISmith & Contill2QQ'8[) powerful winds of the high-mass stars afreet the clouds m 

^ The CNC extends over at least ~ 80 pc, which corresponds to ^^''^ '^'^^'^''^ ^7^-. *e one hand, the process of photoevap- 
^ more than 2 degrees on the sky, demonstrating that wide-field ^^^"§1^ irradiated cloud surfaces can disperse even 

- - surveys are necessary in order to obtain a comprehensive rather massive clouds on relatively short timescales; transform- 
information on the full star forming complex. Several wide-field ^ense molecu ar clouds into warm/hot low-density atomic 

surveys of the CNC have recently been cairied out at different Sf ^ ^"fy ^"^'} *e potential for further star formation, 

wavelengths. The combination of a large Chandra X-ray survey other hand the compression of clouds by irradiation 

(see Townslev et al. 2011b) with a deep near-infrared survey and by expandin g>bbles driven by evolving HII regions (e.g. 

[Preibisch et al. (201 Ic b) and Spitzer mid-infrared observations Deharveng et al. 2005) or stellar winds can produce gravitation- 

dSmith et al. 2010; Povich et al...20in provided comprehensive ^l^^ unstable density peaks and lead to triggered star formation 
information about the young stellar populations. Our sub- Smith^etdJI^lO). The detailed balance of these two op- 

posmg processes determines the evolution of the complex and 
* The //.r.c/../ data described in this paper have been obtained in the ^^^^^^f ^^^^ of the original cloud mass is actually trans- 
open time project OTl_tpreibis_l (PI: T. Preibisch). Herschel is an ^^med into stars and what fraction is dispersed by the feedback. 
ES A space observatory with science instruments provided by European- In order to study this interaction between the stars and the sur- 
led Principal Investigator consortia and with important participation rounding clouds, a comprehensive characterization of the cloud 
from NASA. structure and temperature is a fundamental requirement. 
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The ESA Herschel Space Observatory (iPilbratt et al. 11201 Oh 
is ideally suited to map the far-infrared emission of the warm 
and cool/cold molecular clouds and is currently observing many 
galactic (and extragalactic) star forming regions. We have used 
the Herschel Observatory to perform a wide-field (^5 square- 
degrees) survey of the CNC (RL: Th. Preibisch) that covers 
the entire spatial extent of the clouds. Far-infrared photomet- 
ric maps were obtained with PACS (Poglitsch et al. 2010) at 70 
and 160 yum, and with SPIRE ( Gri ffin et al. 1 1201 Ol) at 250, 350, 
and 500 yum. A first analysis of these Herschel observations, fo- 
cusing on the larg e scale structu r es and global properties, has 
been presented in Preibisch et al. (2012) (Paper I in the follow- 
ing text). A detailed study of the young stellar and protostel- 
lar population detected as point- s ource s in the Herschel maps 
is described in iGaczkowski et aP (l2013l) . The properties of the 
clouds around the prominent HII region Gum 31, located just 
north-west of the central Carina Nebula, as derived fr om these 
Hersch el and other observations are studied in O hlendorf et al.l 
(l20T3h . In this paper, we present a detailed analysis of the tem- 
perature and column density of the small-scale structures of the 
entire CNC. 

This paper is organized as follows: in Section |2] we briefly 
summarize the Herschel PACS and SPIRE observations, data re- 
duction and the final temperature and column density maps of 
the CNC. In Section[3]we compute the temperature, column den- 
sity and UV-flux maps of the clouds. In Section |4l we present a 
detailed analysis of the cloud structures resolved by Herschel in 
the CNC, e.g. the southern pillars. We discuss, in Section |6l the 
comparison between dense and diffuse gas, the relation between 
dust and gas mass estimates, the strength of the radiative feed- 
back and the CNC as a link between local and extragalactic star 
formation. A summary of our results and conclusions is given in 
Section [TJ 

2. The Herschel PACS and SPIRE maps of the 
Carina Nebula Complex 

2.1. Observations and Data Reduction 

The CNC has been observed on December 26, 2010 using the 
parallel fast scan mode at 60'7s, obtaining simultaneously 70 fim 
and 160 fim images with PACS, and 250, 350, and 500 yum maps 
with SPIRE. 

The scan maps cover an area of 2.8° x 2.8°, which corre- 
sponds to about 1 10 pc X 1 10 pc at the distance of the CNC. The 
data reduction has been carried out combining the HIP E (v. 7.0; 
Iptt 2010) with the SCANAMORPHOS package (v. 13.0: iRoussel 
|2012). All the details of the reduction can be found in Paper 
I. We here only emphasize that we used the option galactic 
in SCANAMORPHOS which preserves the brightness gradient over 
the field. The final angular resolution of the FACS/Herschel 
maps was ^ VI" - W\ while for the SFIRE/ Herschel maps 
^ 20" - 36". 

2.2. Calibration 

The observations and mosaics of the CNC have been previously 
presented in Paper I. In this work an additional eff'ort to cali- 
brate the final mosaic has been done by discussing the eff'ect of 
the large scale flux loss and taking into account the color correc- 
tion, which might influence our further analysis. In the first case, 
it is important to recall that Herschel is a warm instrument, not 
designed to obtain absolute fluxes. The warm emission from the 
telescope's mirror is removed using a high-pass filtering during 



the data reduction. However, this step might remove some large 
scale emission in the final mosaics. The flux loss can only be es- 
timated by extrapolating the Herschel fluxes to the wavelengths 
observed by a cold telescope, such as e.g. Planck. 

Another possibility is to obtain the contributions from back- 
ground and foreground material estimating the dust in fore- 
ground/background gas that is traced by H I, for the atomic emis- 
sion, and CO, for the molecular emission. Th is approach has 
been presented by iRivera-In graham et al.l (l2013b . 

Since the Planck data are not yet publicly available and we 
do not have H I and CO spectral data of the CNC, we compared 
our Herschel 70 jim map to the IRAS 60 jim map. This approach 
was already used in Paper I and the details can be found there. 

We also use another approach, which consists in removing 
the local background at each wavelength, by assuming an aver- 
age value over the field (see e.g. Stutz et al. 2010). In this way 
we only take into account the Herschel emission above this local 
background. This was computed in a region of the PACS and 
SPIRE maps outside the CNC where there was no cloud emis- 
sion. We find that the local background is 0.0005 Jy/arcsec^. 

We checked whether the subtraction of the local background 
estimated was aff'ecting our result on the temperature structure. 

Both approaches show that the large scale flux loss can intro- 
duce an error negligible compared to a conservative calibration 
error of 20% of the flux adopted in our work. 

The second calibration step taken into account is the color 
correction. This correction can afl'ect the Herschel photometry 
at the shortest wavelengths in the case of low temperature^jTo 
correct the PACS fluxes we used the first temperature map ob- 
tained by fitting the SED on the uncorrected Herschel fluxes 
(Func)- For a given temperature, we applied the color corrections 
ccpacsio and ccpacsieo listed in the PACS manual^ the corrected 

flux Fcorr is COmputcd as Fcorr = FunclcCpacs- 

In a second step the fit of the SED has been repeated using 
the color-corrected fluxes. 

Between the two iterations the temperature map difl'ers of 
about 0. 1 K in the coldest part of the CNC (T~ 16 K). For warmer 
parts the correction was less than 0.1 K. 

We checked the efl'ect of the color-correction on the SPIRE 
photometry. The SPIRE color corrections are obtained checking 
the slope a of the spectral energy distribution (SED hereafter) 
between 250 and 500 fim. From an average spectral index value 
of about 2.4 over the entire Herschel mosaic, the color correction 
would have ranged between 1.02 and 0.970. This correction is 
hence negligible. 

3. Analysis 

3.1. Temperature and column density of the clouds 

In this section we describe in detail the procedure to obtain a 
precise temperature and column density map of the CNC. The 
most reliable and powerful approach to derive the dust temper- 
ature and density maps is fitting, pixel by pixel, the SED to the 
fluxes in the 5 Herschel bands. The final Herschel mosaics at 
70, 160, 250, 350, and 500 yum (derived as previously described) 
have been convolved to the angular res olution of the 500 jim im- 
age, using the proc edure p resented in Stutz et al.l (1201 Oh and the 
kernels from iAniano et al.l (1201 ID . 

^ PACS manual: "PACS Photometer Passbands and Colour Correc- 
tion Factors for Various Source SEDs" 

^ SPIRE color corrections from the SPIRE Photometry Cookbook; 
Bendo et al. (2011) 
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Table 1. Values of the dust mass absorption coefficient, ky interpolated 
at the Heschel bands from the model of Ossenkopf & Henning, (.1994) 
with a p value of 1.9. 



A 70yum 


160yum 


250 yum 


350yum 


500 yum 


^^[cmVg] 118.0 


24.8 


11.6 


5.9 


2.9 



Sinc e the nebula is opticall y thin at all the Herschel wave- 
lengths (iPreibisch et al.ll2011bl) . the emission of the nebula is 
By(T) ' T, where By is the Planck function at a temperature T 
and T is the optical depth, which is proportional to the dust 
mass absorption coeffi c ient, k y. From the models presented in 
lOssenkopf & HenningI (Il994l) . we choose a dust model with a 
highest JS value (i.e. 1.9) which is representative to describe 
a dense molecular cloud. This model has a standard MRN- 
distribution for diffuse interstellar medium ( Draine & Leell 19841) 
composed of grains without ice mantels. We used the tabulated 
values of ky [cm^g] for such a model and interpolated at the 
Herschel wavelengths. In Table [T] we show the values of ky that 
we obtained. 

The fit of the SED is obtained by using a black-body leav- 
ing as free parameters the temperature T and the surface density 
2 [g/cm^]. A similar approach to obtain temperature and col- 
umn density has been applied f or the anal ysis of the HOBYS 
key project by e.g. Hennemann et alJ (l2012h and the Gould Belt 
key program by e.g. Palmei rim et al.l (l2013l) . 

The final temperature and density maps are obtained from 
the Herschel/FACS and SPIRE color corrected mosaics. 

The final temperature map of the CNC is shown in Figure [T] 
We find that the average temperature of most of the nebula is 
about 30 K, ranging between 35-40 K in the central clouds and 
26 K in the clouds at the edge of the nebula. 

In Paper I we presented a first color-temperature map of the 
CNC, derived from the ratio of the 70yum and 160yum images. 
The warmer cloud surface temperature is traced with the best 
resolution, since the convolution of the two original mosaics has 
been done on the PACS 160 jim. 

However, a color temperature may be less well suited for 
dense clouds where much of the 70 yum - 160 yum emission comes 
from the warmer cloud surfaces. In this case, the temperature 
computed by the SED fitting computes the beam averaged tem- 
perature among the line of sight and it is more sensitive to the 
densest (and thus coolest) central parts of clouds. 

We find that the temperature computed in Paper I show val- 
ues up to ~ 5% higher than the values with the SED fitting. This 
result was already expected because the color temperature is bi- 
ased to the warmer cloud surface. 

The second free parameter of the SED fitting is the dust sur- 
face density, X. The column density Nr is obtained as following 

A^H = 2Nn, = ^'^'^ (1) 

where mn is the hydrogen mass and jiu2 is the mean molecular 
weight (i.e. 2.8). Multiplying by the gas-to-dust mass ratio R 
(assumed to be 100), we obtain the total column density. 

In Figure |2] we show the resulting total column density 
map. The median value over the entire CNC (2.3°x2.3°) is 
1.8 X 10^^ cm"^, while the mean value is 2.3 x 10^^ cm"^. 

We find column densities up to 5 x 10^^ cm"^. From the 
Southern Pillars, in the South-East direction from Carina, there 
is an almost continuos structure with density of about 3.0 - 



5.5 X 10^^ cm ^. In the opposite direction, north-east to south- 
west from the Carina Nebula, the density is lower, down to 
9x 1020cm-2. 

In particular the large elongated bubble south-west of r] Car 
is a quite prominent empty region already described in Paper I. 

We note that a conservative calibration error of 20% in the 
intensity causes just 5-10% uncertainty in the derived temper- 
atures, and at T~20 K, a 5% error in temperature will cause a 
<10% error in the derived column density. 

3.2. Determination of the UV-flux 

The surfaces of the clouds we see in the FIR maps are mainly 
heated by the far-UV (FUV) radiation from the hot stars in the 
complex; the FUV photons are very efficiently absorbed by the 
dust grains and thus heat the dust to temperatures well above 
the "general ISM" cloud t emperatures of ~ 18 K (see, e.g. 
iHoUenbach & TielenslfT999h . Since the heated dust grains cool 
by emitting the FIR radiation we can observe with Herschel, one 
can estimate the strength of the FUV radiation field from the 
observed intensity of the FIR radiation. This yields important 
quantitative information about the local strength of the radiative 
feedback from the massive stars. 

For a few selected clouds in the CNC, the strength of the 
FU V irradiation has bee n already determined in previous stud- 
ies. iBrooks et"aD (l2003l) investigated the photo-dissociation re- 
gion (PDR) at the eastern front of the dense cloud west of Tr 14 
and found that the FUV field is 600 to 10 000 times stronger than 
the so-called Habing field, i.e. the average intensity of the local 
Galactic FUV ( 912A-2400 A) flux of 1.6 x 10"^ ergcm'^ s"^ 
(iHabingll 19681) . iKramer et al.l (i2008) derived an independent es- 
timate of the FUV irradiation at the same cloud front, finding 
values of 3400 to 8500 times the Habing field. For for another 
cloud south of T] Car, these authors found irradiation values of 
680 to 1360 times the Habing field. 

Us ing the approach described in more detail by lKramer et al.l 
(l2008l) , we computed the FUV radiation field for all clouds in the 
CNC from the observed intensity of the FIR radiation in our Her- 
schel maps. For this, we used our PACS 70 yum and 160 yum band 
maps to determine the total intensity integrated over the 60 fim to 
200 yum FIR range. Since these two bands cover the peak of the 
FIR emission spectrum of the irradiated clouds, they constitute a 
good tracer of the radiation from the heated cloud surfaces, and 
also provide very good angular resolution. We did not include 
the longer wavelength SPIRE bands, firstly because the angular 
resolution in these bands is considerably worse, and secondly 
because the emission at ^ 250 yum is probably dominated by the 
thermal emission from cold dense clumps, that are well shielded 
from the ambient FUV field; the emission from these very dense 
and cold cloud structures does thus not directly trace the FUV 
irradiation. In any case, the contribution of the SPIRE bands to 
the integrated intensity would be quite small. 

To determine the total 60 - 200 yum FIR intensity, /fir, we 
first smoothed the PACS 70 yum map to the angular resolution of 
the PACS 160 jim map and created maps with a pixel size of 4.5'' 
(corresponding to 0.049 pc). For each pixel in this map, we then 
multiplied the specific intensity in the two bands by the band- 
width (60 - 80 yum and 125 - 200 yum; see PACS MANUAlB- In 
order to fill the gap from 80 yum to 125 yum, we took the average 
of the 70 yum and 160 yum intensities and multiplied it by the cor- 
responding band width. The resulting sum is a measure of the 
total 60 - 200 yum FIR intensity (/fir). 

^ http ://herschel. esac . esa.int/Docs/PACS/html/pacs_om.html 
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Fig. 1. Temperature map of the Carina nebula. The dashed boxes represent: the central region around t] Car (box 'a'), the Southern Pillars region 
(box 'b'), the Treasure Chest Cluster (box 'c'), and the 'wave' pattern (box 'd'). These regions are analyzed in detail in Figures |4][5l 13 and[8] The 
center of the figure corresponds to the position {anooo, ^moo) - (10^ 45"^ 2P, -59° 34' 20"). 



Following the arguments of lKrameretaP (l2008h . the FUV 
field strength can then be computed via the relation 



Go 



47r/FiR 



1.6 X 10"^ erg cm"^ s" 



(2) 



where the denominator contains the Habing field. 

The resulting map of the derived FUV field in units of the 
Habing field is shown in Figure [3] In the enlargements of the 
FUV map, shown in Figures IH [S] and [71 we also overplotted 
the positions of the known high-mass stars in the CNC. Fo r the 
individual clouds already studied by iKramer et al. 1 (HOOS) and 
iBrooks et al. (2003), our values agree well with the previous, 
independent determinations by these authors. 

Our map shows that nearly all cloud surfaces in the CNC 
are irradiated with Go values above ~ 300, which seems rea- 
sonable f or a regio n containing such a large number of OB stars 
(j Hollenbach et al.Ll991) . Particularly high Go values, between 
^ 3000 and ^ 10000, are found for the massive cloud to the 
west of Tr 14, Values above 3000 are found in the Keyhole Neb- 
ula, several clouds to the north and west of Tr 14, and a pillar 

south-east of ri Car. 

It is interesting to compare our map to the results of ISmithI 

(l2006l) . who computed an estimate of the radial dependence of 
the average FUV radiation field from the spectral type of each 
high mass star in the central clusters (i.e. Tr 14, Tr 15 and Tr 16). 
At distances larger than 10 pc from the center of the Nebula, the 
ionizing source can be treated as a point-like source, and the 



ionizing flux drops as r"^. The local radiation within few par- 
sec from each high mass cluster is found to be more important 
than the cumulative effect from all the clusters. The FUV radia- 
tion computed by Smith ( 2006) is consistent with the values we 
found: between 100 and 1000 at distances larger than 10 pc from 
the center, and values up to 10"^ -10^ near the individual clusters. 

In the Southern Pillars area, which is further away from the 
hot stars in the central clusters, we find Go values around 1000 
and up to 2000 along the pillar surfaces. 

In most parts of our map area, the spatial distribution of the 
derived FUV intensities correlates well with the cloud tempera- 
tures. However, in some locations, clear diflferences can be seen. 
One prominent case is the Gum 3 1 region in the north-western 
part of our maps (see lOhlendorf et all2013i for a comprehensive 
study of this region). While the peak of the cloud temperature is 
found in the difl'use gas around the OB stars in the central cluster 
NGC 3324, the highest FUV intensities are found at the western 
and southern rim of the bubble surrounding NGC 3324. 



4. Detailed structure of the clouds in the Carina 
Nebula Complex 

In this section we characterize particularly interesting small- 
scale structures in the CNC which are resolved in the Herschel 
maps. In particular we will analyze the central part of the CNC, 
the Southern Pillars region and the Treasure Chest Cluster . 
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Fig. 2. Column density map of the Carina nebula. The color scale is expressed in cm~^. The dashed boxes represent: the central region around 
T] Car (box 'a'), the Southern Pillars region (box 'b'), the Treasure Chest Cluster (box 'c'), and the 'wave' pattern (box 'd'). These regions are 
analyzed in detail in Figures |4l|5ll2l and[8] The center of the figure corresponds to the position (aj2ooo» ^moo) = (10^ 45"^ 21% -59° 34' 20''). 



4.1. Central part 

An enlargement of a 23 pc x 26 pc region around rj Car is shown 
in Figure IH Our Herschel maps provide sufficient angular res- 
olution to investigate the temperature and column density struc- 
ture of the clouds around r] Car in detail. 

In the southern direction from ri Car a dense molecular cloud 
ridge is resolved in the IRAC image and column density map. 
The temperature is not homogeneous around t] Car: in the North- 
West direction from ri Car the temperature is higher, ranging be- 
tween 35 and 45 K; in the South-East direction the temperature 
is lower, between 25 and 35 K. 

In this region, there are the young clusters, Trumpler 14, 15 
and 16, which host about 80% of the high mass stars of the entire 
complex. This is also the hottest region of the nebula with tem- 
peratures ranging between 30 and 50 K (excluding the position 
of ri Car itself which reaches values of 60 K). 

In particular, the clouds in and around the stellar cluster 
Tr 16, which hosts most of the high mass stars in the nebula 
including r] Car, have temperatures between 29 and 48 K. Tr 14, 
which hosts 20 0-type stars, shows cloud temperatures of 30- 
33 K, and Tr 15, with 6 0-type stars, displays a lower temper- 
ature of 26-28 K. The local temperature of the cloud is hence 
related to the number of high mass stars. 

The column density map of the molecular cloud around 
the clusters shows low values of about 9.0 x 10^^ cm"^. The 
molecular cloud structures at the edges of the clusters are colder 
(~ 23 K) and denser ~ 3.9 x 10^^ cm"^ compared to the central 



part. The most prominent cloud West of Tr 14 has a temperature 
of about 30 K and a decrease in density from the inner to the 
edge part. 

4.2. The Southern Pillars 

In the Southern Pillars region of the CNC, the temperature 
ranges between 21 K and 30 K. For a few of the most prominent 
pillars, we analyze the temperature and density profiles in detail. 
In Figure[5] an image of the Southern Pillars is compared with 
the temperature and column density maps. All the structures in 
the composite image in Figure [5] are reproduced in the density 
map. The two most prominent pillars (Pillar A and Pillar B in 
Fig. O have a projected angular separation of about 4'\ In the 
East direction from Pillar A and Pillar B another pillar is evident 
in the composite image of Figure [51 This is only barely visible 
in the density map, while in the temperature map only a slightly 
colder temperature compared to the local environment is seen at 
the head of the pillar. The projected size of both. Pillar A and 
Pillar B, is about 5 pc. The head of the pillars points in the direc- 
tion of the central stellar cluster Tr 16. Their column densities 
range from ~ 3.0 x 10^^ cm"^ at the edge of the pillar structure, 
up to ~ 10x10^^ cm"^ at the peak of density in the central part of 
Pillar B. In order to analyze the temperature and density profiles 
along the pillars, we show the run of temperature and column 
density along cuts perpendicularly to their main axis, spaced by 
~r from the bottom to their head. While the temperature pro- 
file always decreases from the edge to the center of the pillar. 
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the density profile increases and the maximum of the density in 
the pillar's center corresponds to the minimum temperature. The 
temperature varies by 2 K according to their position on the pil- 
lars: cooler at the bottom of the pillar and hotter at their head. 
The Pillar A is found to be about 2 K cooler than Pillar B. A 
possible explanation is that Pillar B is closer to the stellar cluster 
Bochum 1 1 which harbours 5 high-mass stars. The temperature 
and density profile of Pillar A suggests that it consists of two 
partly resolved pillars: the first and third sections show a double 
peak in both temperature and density profiles, while the second 
and fourth sections show only one peak which traces one of the 
two here unresolved pillars. 

Integrating our column density map over the spatial extent of 
the entire "giant pillar" (i.e. including both, Pillar A and Pillar 
B), we derive a total mass of ^ 9350 M© for the pillar structure. 

The density structure of the p illars ha s been also pre- 
dicted by theoretical s tudies, e.g. by iGritsctm eder et al. (2010); 
iTremblin et al.l (|2012h . In the simulations of Grit schneder et all 
(12010) Mach 43, after 500 kyr the ionisation leads the cre- 
ation of few big p i llars. The most prominent pillar created by 
iGritschneder et al.l (12010*) has a size o f 0.4 pc x 0.8 pc , with col- 
umn density up to 1.3 x 10^^ cm'^. In lTremblin et al.l (f2Q12) the 
formation of a slightly bigger 1 parsec long pillar takes place in 
the Mach 1 simulation between 200 kyr and 700 kyr. In partic- 
ular their non-turbulent simulation was based on the high curva- 
ture curvature of the dense shell which leads to pillars. 

The size s of the pillars p redicted by both, Gritschneder et al.l 
(120 Id) and iTremblin et all (l2012h . are between four and two 
times smaller compared t o the extensions of P illar A and Pil- 
lar B. In the simulation of ITremblin eTaP (l2012h It IS mterestmg 



^ Details of the Mach 5 simulations are in IGritschneder et aTl (l2009h . 



to notice that the pillar has two dense heads separated by 0.2 pc. 
This might explain the double peak in the temperature and den- 
sity profile found in Pilla r A. The column density predicted by 
IGritschneder et al.l (1201 Ol) is in good agreement with the profile 
of Pillar B. 

4.3. The cloud around the Treasure Chest Cluster 

In the Southern Pillars region, the Treasure Chest Cluster is 
the largest and most prominent embedded star cluster. It is 
partly embedded in a dense cloud structure and is thought to 
be very young, with an estimated aged of less than 1 Myr (e.g. 
[Smith etal.1 120051: iPreibisch et al.l 1201 Id). The brightest and 
most massive cluster member is the star CPD -59° 2661 (= MJ 
632), for which a spectral type of BlVe is Hsted in the SIMBAD 
database. In the near-infrared image (Figure |7), the cluster is 
surrounded by an arc of diff'use nebulosity, which seems to rep- 
resent the edge of a small cavity that opens to the south-west. 
About 30'' north-east of the Treasure Chest Cluster, and clearly 
outside this cavity, another bright star. Hen 3-485 (= MJ 640, 
spectral type Bep) is seen. 

Our temperature, column density, and FUV intensity maps 
of the cloud around the cluster show a very interesting morphol- 
ogy, even if the structures are not well resolved. The compar- 
ison of the near-infrared image with our Herschel maps shows 
that the Treasure Chest Cluster is situated near the head of a 
pillar-like cloud. Our column density map shows that the clus- 
ter is clearly off'set from the density peak (2.3 x 10^^ cm"^) of 
the pillar. The map also reveals a "kidney"-shaped column den- 
sity depression near the location of the cluster. Our temperature 
map shows colder gas (7 25 K) near the column density peak, 
and considerably warmer gas (T ^ 32 K) at the position and to 
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Fig. 4. Upper left: Composite optical DSS (in blue), Spitzer IRAC4 (in green), and LABOCA (in red) image of the central part of the Carina 
Nebula Complex around r] Car (box 'a' in Fig.[T]and|2]). Upper right: Temperature map of the central part of the Carina Nebula Complex. Lower: 
Column density map, expressed in cm~^ {left) and Far-ultraviolet (FUV) flux Go (in units of the Habing field) {right) of the enlargement around 
77 Car. The blue crosses represent the positions of the high mass stars in this region and the black star the position of 77 Car. 



the south-west of the cluster, as expected due to the local heat- 
ing of the dust and gas in the cavity by the cluster stars. These 
results support our interpretation that the cluster is located in a 
small cavity at the western edge of the pillar, which is open in 
the south-western direction. 

The star Hen 3-485 is seen (in projection) almost exactly 
at the column density peak of the cloud. However, its optical 
brightness and the moderate redenning (the observed color of 
B - y = 0.7 suggest a visual extinction of Ay ^ 3 mag) show 
that it is most likely located in front of the pillar. 

Our map of the FUV intensity shows a clear peak just be- 
tween the two B- stars, coinciding well with the arc of reflection 
nebulosity in the near-infrared image. 

From our column density map we derive a total mass of ^ 
780 Mq for the pillar-like cloud associated to the Treasure Chest 
Cluster. 

5. Detection of a "wave"-like pattern 

In the region between the central parts of the Carina Nebula and 
the Gum 31 region, in the northern part of our Herschel maps, 
our column density map reveals an interesting pattern of quite 
regularly spaced "waves" (see also Figure 2 and 3 of Paper I). It 



consists of a pattern of parallel lines (see Figure [S} running to- 
wards the north-eastern direction. The wavelength of this pattern 
is 14.7', corresponding to 9.8 pc. The column density varies 
from ~ 1.8 X \(fi^ cm"^ at the maxima to < 9.0 x 10^^ cm"^ at 
the minima. 

We first checked whether such a "wave"-like pattern can be 
an eflTect of the observation mode of Herschel. The observa- 
tions are performed with two scanning directions which are later 
combined into the final map. The first scan direction is at 8° and 
the second one is at 100°, both from West to North direction. 
The "wave"-like pattern is found instead at 40° from the West to 
North direction, which does not correspond to the scan direction. 
We also created a final image from each scan direction and the 
"wave"-like pattern is present in both of them. For this reason 
we exclude that this is an observational eflTect. 

One possible explanation for this pattern is that we see waves 
at the surface of the molecular cloud, similar to the pattern dis- 
covered in the Orion Nebula by IB erne et al.l (1201 Oh . 

In the case of the Orion Nebula, the waves have been ex- 
plained by a Kelvin-Helmholtz instability that arises during the 
expansion of the nebula as the gas heated and ionized by mas- 
sive stars is blown over pre-existing molecular clouds that are 
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Fig. 5. Upper left: Composite image (LAB OCA 870 yum, Spitzi 
Fig. [T] and 12]). Temperature map {upper right). Column density m 
of the Habing field) {lower right) of the enlargement of the regie 
temperature and column density are shown in Figure |6] 




Delta R.A. [arcmin] 

• 3.6 and 8.0 jjm) of the giant pillar in the Southern Pillars region (box 'b' in 
p, expressed in cm~^ {lower left) and Far-ultraviolet (FUV) flux Go (in units 
[. The solid lines highlight the cut lines of the prominent pillars for which 
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Fig. 7. Upper panel: Composite image of the Treasure Chest Cluster (box 'c' in Fig.[T]and[2l): HAWK-I /-band image in blue, HAWK-I ^^-band 
image in green, IRAC3 image in red, and the (red) contours of the LAB OCA map. In the following panels are shown the temperature {upper 
right), column density maps {lower left) and Far-ultraviolet (FUV) flux Go {lower right). The column density map is expressed in cm~^, while the 
FUV in units of the Habing field. The black 'X' represents the position of the Treasure Chest Cluster, while the black cross the position of the star 
Hen 3-485. 



located a few parsecs away from the center of the cluster, where 
the young massive stars are concentrated. 

The situation in the Carina Nebula is similar: The feedback 
from the numerous massive stars in t he central regions of the 
Carina Nebula drives a bipolar bubble (ISmith et al.ll2QQQh . 

The location of the observed wave pattern agrees roughly 
with the northern part of this bubble system. The waves are 
seen at a projected distance of ~ 6 - 15 pc from the massive 
stars in Tr 16 and Tr 14. For a comparison to the Orion Neb- 
ula case, this larger distance may be easily compensated by the 
much larger number of massive stars in the Carina Nebula and 
the correspondingly higher level of stellar feedback. Therefore, 
the density and velocity of the hot plasma that is thought to cause 
the wave pattern in the Orion Nebula may well be of a similar 
order of magnitude as in the pattern we see in the Carina Nebula. 

Some information about the outflowing material can be ob- 
tained from the Chandra X-ray survey in the conte xt of the 
Chand ra Carina Complex Project (CCCP) (see Townsley et al.l 
l2Qllbh . A detailed analysis of th e diflFuse X-ray emission , which 
traces the hot plasma, is given in lTownsley et al.l (l2Qllah . 



The largest part of the wave pattern revealed by Herschel lies 
just outside the area covered by the Chandra X-ray map. At the 
north-eastern edge of the Chandra X-ray map, which touches the 
south-western edge of the wave pattern revealed by Herschel, the 
dominant plasma components have temperatures ranging from 
0.14-0.17 keV to 0.25-0.27 keV. 

The ROSAT image, that covers a larger field-of-view than 
the Chandra map, shows that the amplitude of the difl'use X-ray 
emission in the area of the wave pattern is rather weak. This 
suggests that the hot plasma in this region is not confined, but 
streaming outwards. 

The most remarkable difiTerence to the waves in the Orion 
Nebula is the substantially longer wavelength of the wave pattern 
in the Carina N ebula, ~ 10 pc rath er than ~ 0.1 pc. According 
to the models of lBeme et al.l (12010 ), this longer wavelength sug- 
gests that a plasma of density ~ 1 cm"^ with a flow velocity of 
-100 km/s could be causing the instability. Such values appear 
quite possible in the case of the Carina Nebula superbubble. 
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Fig. 8. Temperature (left) and column density, expressed in cm ^, (right) map of the north-eastern part of the Carina Nebula Complex (box 'd' 
in Fig. [Hand 12]). The dashed lines represent a "wave pattern" which seems to connect the two main regions of the complex. 



6. Discussion 

6.1. Comparison of dense and diffuse gas 

In Fig. [9] we compare the Hersc hel 160/immapto the LAB OCA 
870 yum map presented in Preibi sch et al.l (1201 Idl) . This combi- 
nation has the advantage that the angular resolution of the Her- 
schel 160 yum map (^ 12" - W) matches the angular resolution 
of the LABOCA map (18'0 quite well. The differences in these 
maps are thus not due to the difference in angular resolution (as 
would be the case for the other Herschel maps) but rather trace 
the differences in the spatial distributions of the very cold gas 
and the somewhat warmer gas. 

In most clouds, the cold gas (traced by the sub-mm emission) 
is considerably more spatially concentrated than the warmer gas 
(traced in the Herschel 160yum map). The most prominent ex- 
ample is the dense cloud west of Tr 14, in which the cold dense 
cloud center is surrounded by a much more extended halo of 
widespread warmer gas. 

However, in many of the pillars we find a very different sit- 
uation. Most pillars that are seen ~ 15' - 20' to the south-east 
of ri Car, are not surrounded by significant amounts of extended 
warm gas. The spatial extent of the warmer gas matches that of 
the cold gas quite well. 

The pillars at the southern periphery (more than ~ 20' south 
of 7/ Car) show filaments of warmer gas at their base that seems to 
stream away from the densest parts of these clouds in the direc- 
tion opposite to the pillar head, just as expected in the scenario 
of cloud dispersal by ionizing irradiation. 

A particularly interesting difference in the spatial configura- 
tion of the cold gas versus the warmer gas is displayed by the 
cloud structure ~ 5' - 10' to the south and south-east of ri Car 
(that constitutes the "left wing" of the V-shaped dark cloud). 
Here, the cold gas is highly concentrated in two distinct com- 
pact clouds (labeled as 'A' and 'B' in Fig.O, which have a gap 
(with a width of about 3') between them, just at the position 
that is closest (in projection) to ri Car. In a remarkable contrast, 
the configuration of the warmer gas resembles just a continuous 
elongated distribution, with no indication of any gap at this posi- 
tion. This implies that there is a considerable amount of warmer 
gas in the region between these two clumps of cold gas. 



6.2. Dust vs. molecular gas mass estimates for individual 
clumps 

The area covered by our Herschel maps includes 14 of the C^^Q 
clumps detected and characterized in the study of lYonekura et al.l 
( 2005). Compact clouds are seen in our Herschel maps at the 
locations of all of these clumps with exception of C^^O clump 
number 10. 

In order to derive clump mass estimates from their FIR emis- 
sion, we integrated over boxes centered at the cloud positions 
as seen in our Herschel maps th at correspond to the re ported 
positions of the C^^O clumps of lYonekura et al.l (l2005b . Our 
box size s were chosen to be a s similar as possible to the sizes 
given bv lYonekura et al.1 (l2005l) (details in Table[2]). Eight of the 
C^^O clumps are also located in the field of view of our previous 
LABOCA observations (Preibisch et al. 2011c), from which we 
already estimated clump masses based on their sub-mm fluxes. 
In TableOwe summarize the total, gas + dust, clump masses and 
temperatures computed using the different methods. 

This comparison shows that the masses computed from the 
our Herschel maps agree generally quite well to the previous es- 
timates based on the sub-mm fluxes. However, the masses based 
on FIR and sub-mm data are alwa ys about 2-4 times lower than 
the masses estimated by Yonekura et al.l (l2005l) from their C^^O 
data. 

There are several possible explanations for this discrepancy. 
One major issue are the extrapolation factors used relating the 
measured dust or CO molecule emission to the total gas + dust 
mass. It is well known that the gas-to-dust ratio, for which we 
assumed the canonical value of 100, is uncertain by (at least) 
about a factor of 2. Concerning the masses derived from CO, it 
is known that the fractional abundance of carbon-monoxide may 
not be constant but depend on the cloud density; thus, the factor 
used to extrapolate the CO mass to the to tal mass can easily var y 
by factors of 2-3 (see, e.g., discussion in [Goldsmith et al .120081) . 

We also find small, but systematic differences in the derived 
cloud temperatures. In most clumps, the dust temperature we de- 
rive from our Herschel maps is a few Kelvin higher than the CO 
gas temperatures derived by |2] This is probably related to the 
fact to both temperatures represent averages over (i) the corre- 
sponding beam area and (ii) along the entire line-of- sight though 
the clumps. Since real clumps are not spatially isothermal but 
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Fig. 9. Herschel 160yL/m map (in Jy/px) with the LABOCA 870yL/m contours overlaid at 0.08, 0.16, 0.32, 0.64, 1.28, 2.56 Jy/beam. Highlighted 
is the position of 77 Car. 



have some kind of (unresolved) density and temperature struc- 
ture, the different biases of the different observations result in 
these different temperature values. 

We finally note that the Herschel maps show many more 
dense cold dusty clouds, which remained undetected in the C^^O 
maps. 

6.3. Cloud mass budgets based on different tracers 

In Table [3] we summarize the integrated cloud masses for differ- 
ent parts of the CNC derived from different tracers. The masses 
derived in this paper via fitting of the FIR spectral energy distri- 
bution, and in iPreibisch et al.l (1201 id) from sub -mm fluxes, are 
based on dust as a mass tracer. The study of lYonekura et al. l 
(|2005) gives mass estimates derived from different CO isotopo- 
logues, and is thus based on gas tracers. 

Considering the total cloud complex, we find rather good 
agreement between the different mass estimates for specific col- 
umn density thresholds. Our Herschel based cloud mass esti- 
mate for an column density threshold corresponding to Ay > 1 
is well consistent with the mass estimate based on ^^CO. This 



is in good agreement with the standard idea th at molecules ex- 
ist in side the shielded regions of dark clouds dGoldsmith et al.l 
[2QQ8h. The threshol ds in Ay of 0.9 and 7.3 have been chosen as 
in lLada et al.l (l2012h . 

Our previ ous cloud mass estimate based on the LABOCA 
sub-mm flux (iPreibisch et al. 2011c) agrees well with the mass 
estimate based on ^^CO. This fits with the assumption that the 
LABOCA sub-mm map traces only dense clouds and that ^^CO 
is a tracer of relatively dense gas. 

The total mass of all individual clun ips detected in our 
LABOCA sub-mm map (as determined by iPekruhl et"ani2013h 
is between the cloud mass estimates based on ^-^CO and C^^O, in 
good agreement with the idea that C^^O traces denser gas than 
i^CO. 



6.4. The strength of the radiative feedback 

As mentioned above, the cloud surfaces in the CNC are strongly 
irradiated: we find Go values generally between ~ 1000 and up 
to ^ 10000. These numbers provide important quantitative in- 
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Table 2. Mass and temperature derived for the C^^O clumps from lYonekura et alldlOOSh . 
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References. This work; ^'^^ iPreibisch et al.l(l2011ch assuming a temperature of 20 K; ^^ Wonekura et al.l(l2005h . scaled to the distance of 2.3 kpc. 



formation on the level of the radiative feedback from the massive 
stars on the clouds. 

To put these numb ers in context, we refer to the compi lations 
in lBrooksetaP (l2003h (Tab. 3) and Schneider & Brooksl (120041) 
(Tab. 2), where values of the FUV field strengths in photo- 
dissociation regions (PDRs) in diff'erent star forming regions are 
listed. This comparison shows that the FUV field strength at the 
cloud surfaces in the CNC is very similar to that measured for 
PDRs in the 30 Dor cluster, which is generally considered as the 
most nearby extragalactic starburst system. This result supports 
the idea that the level of massive star feedback in the CNC is 
similar to that found in starburst systems. 

The comparison also shows that, sometimes, even higher val- 
ues of the FUV field strength are found at specific locations in 
much less massive star forming regions. For example, in the 
PDR in the Orion KL region. Go values up to 10"^ - 10^ have been 
determined (see references in Schn eider & Brooksll2004l) . Such 
high local values are, however, only found in clouds located very 
close (less than 1 parsec) to individual high-mass stars, i.e. re- 
stricted to very small cloud surface areas. No such extremely 
high irradiation values are seen in the CNC, because the clouds 
are generally at least a few parsecs away from the massive stars. 
This is probably a consequence of the fact that many of the mas- 
sive stars in the CNC are several Myr old and have thus already 
dispersed the clouds in their immediate surroundings (by their 
high level of feedback). 

Despite the fact that some lower-mass star forming regions 
show peaks of FUV irradiation that can locally exceed the lev- 
els we see in the CNC, the total amount of hard radiation (and 
thus the total power available for feedback) is much higher in the 
CNC. 

The high level of the FUV irradiation may also aff'ect the cir- 
cumstellar matter around the young stellar objects in the CNC. 
The irradiation causes increased external heating and photoe- 
vaporation of the circumstellar disks and can disperse the disks 
within a few Myr. This process has been directly observed in the 
so-called ''proply ds" near the most massive stars in the Orion 
Nebula Cluster (O 'dell et al.l[T993l: |Bally et al.J|2003; iRicci et al.l 
[2OO8) and has been theoretically investigated by numerical simu- 
lation s (e.g., iRichling & Yorkell2000l:lAdams et al ]l2QQ4l:IClarkd ' 
I2007h 



Using near-infrared excesses as an indicator of circumstel- 
lar disks, our previous studies of X-ray selected young stars in 
the diff'erent populations in the CNC actually showed that the 
infrared excess fractions for the clusters in the CNC are lower 
than those typical fo r stellar clusters of sim ilar age with lower 
levels of irradiation (iPreibisch et al.ll20TTcl) . This is consistent 
with the idea that the high level of massive star feedback in the 
CNC causes faster disk dispersal. Similar indications have been 
found in other high-mass star forming regions, where the disk 
fract ions were also found to be lower close to high mass stars 
(e.g. lFang etaDl2012l:lRoccatagliata et al.ll201ll) . 



6.5. On the relative importance of massive star feedback and 
primordial turbulence 

A fundamental question of star formation theory is whether the 
structure of the clouds in a star forming region is dominated by 
primordial turbulence or by massive star feedbac k. 

The recent study of ISchneider et al.l (1201 Oh analyzed Her- 
schel observations of the Rosette molecular clouds and sug- 
gested that "there is no fundamental diff'erence in the density 
structure of low- and high-mass star-forming regions". They ar- 
gue that the primordial turbulent structure built up during the 
formation of the cloud, rather than the feedback from massive 
stars, is determining the course of the star formation process in 
this region. The current star formation process is concentrated 
in the most massive filaments that were created by turbulent gas 
motions. The star-formation can be triggered in the direct in- 
teraction zone between the Hll-region and the molecular cloud 
but probably not deep into the cloud on a size scale of tens of 
parsecs. 

It is interesting to compare the Rosette molecular cloud to 
the CNC. The Rosette cloud has a similar, only slightly lower 
total cloud mass than the CNC, and also harbors a rather large 
population of OB stars. However (as already p ointed out by 
ISchneider & Brooksl 120041: ISchneider etaDll998l) . the Rosette 
Molecular Cloud is exposed to a much (about 10 times) weaker 
UV field than what we find in the Carina Nebula. This order- 
of-magnitude diff'erence in the level of the radiative feedback 
may explain the diff'erences in the cloud structure of these two 
complexes. In strong contrast to the predominantly filamentary 
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cloud structure of the Rosette region, most clouds in the CNC 
show pillar-like shapes. Our Herschel maps reveal (as previ- 
ously suggested on the basis of the Spitzer data; Smith et al. 
I2OIO) the systematic and ordered shape and orientation of these 
numerous pilla rs, that point towards clus ters of high-mass stars. 
As discussed in lGaczkowski et af] (l2Q13l) . the current star forma- 
tion activity in the CNC is largely concentrated to these pillars. 
These results clearly suggest that the cloud structure and the star 
formation activity in the CNC is dominated by stellar feedback, 
rather than random turbulence. 

A strong feedback similar to the situation in the CNC has 
been found in Herschel study of the RCW36 bipolar nebula in 
Vela C bv Minier et al.. (.20131 

6.6. Total FIR luminosity 

From our map of the FIR intensity we also computed the total 
FIR luminosity in the 60 - 200 yum wavelength range, LFiR,tot by 
integrating over a2.3°x2.3° region covering the entire CNC. 
We find a value of LFiR,tot = (1.7 ± 0.3) x 10^ Lq. 

This value is somewhat larger than the total infrared lumi- 
nosity of 1 .2 X 10^ Lq measured bv lSmith & Brooksl (l2007h from 
IRAS and MSX data, but this diff'erence can be explained by 
the sli ghtly smaller area used in the study of ISmith & Brooksl 
(l2QQ7h . that excluded the Gum 31 region (see their Fig. Ic). 

It is interesting to compare this total FIR luminosity to 
the total luminosity of th e known mas sive stellar population 
in the CNC, as derived bv ISmithI (l2006h : Our total FIR lumi- 
nosity is smaller than the total stellar bolometric luminosity of 
2.5 X 10^ Lo, but 70% larger than the total stellar FUV luminos- 
ity of 1.0 X 10^ Lq. The most likely explanation for this apparent 
discrepancy of the total FIR and total FUV luminosity is that 
the current census of massive stars in the CNC is still incom- 
pl ete. Direct su pport for this assumption comes from the study 
of iPovich et al.l (I2OII1) . who could identify 94 new candidate 
OB stars (with Lboi > 10^ L©) among X-ray selected infrared 
sources in the Carina Nebula. They found that the true number 
of OB stars in the CNC is probably about 50%, and perhaps up to 
100%, larger than the currently identified OB population. This 
incompleteness results from the lack of sufficiently deep wide- 
field spectroscopic surveys, that left a significant number of OB 
stars yet unidentified. 

Taking this incompleteness into account, the extrapolated to- 
tal stellar FUV luminosity would match the total FIR luminosity 
we derived quite well. 

6. 7. The Carina Nebula as a link between local and 
extragalactic star formation 

iLada et aP (l2012h collected observational data for numer- 
ous Galactic star forming regions and external galaxies and 
found a well-defined scaling relation between the star formation 
rates and molecular cloud masses above certain column density 
thresholds. They expressed the star formation scaling law for 
these clouds as 

SFR = 4.6 X 10-^ /dg Mg(Mo) M© yr'^ (3) 

where Mq is molecular mass measured at a particular extinction 
threshold and corrected for the presence of helium and /dg is the 
fraction of dense gas. 

They suggest that this fundamental relation holds over a span 
of mass covering nearly nine orders of magnitude. Although the 
exact physical meaning of certain density threshold values is not 



clear (see discussion in iBurkert & Hartmaiinll2012h . this em- 
pirical relation provides an interesting opportunity to compare 
diff'erent star formation regions. In the compilation of obser- 
vational data of ILada et al.l 12012), there is a gap of about four 
orders of magnitude (in cloud mass) between the local Galactic 
clouds and the external galaxies. The CNC is ideally suited to 
bridge this gap and to connect Galactic and extragalactic regions 
in this investigation. 

The cloud masses for the CNC determined with diff'erent 
tracers have been summarized above. For the extinction thresh- 
olds ofAx > 0.1 m ag and Ak > O.H ma^ as used in the study of 
iLada et al.1 (l2012h . we find total cloud masses of ~ 610 000 M© 
and ~ 23 000 M© from our Herschel data of the entire CNC. 

Taking into account the systematic calibration uncertainties 
of the Herschel fluxes as well as uncertainty of the exact bound- 
aries of the CNC, we estimate an uncertainty of 30% for the total 
cloud mass. 

For the star formation rat e in the CNC, several recent deter- 
minations are available. Pov ich et al.l (201 1) analyzed a Spitzer- 
selected sample of YSOs in the central 1.4 square-degrees of 
the Carina Nebula. From this they derived a lower limit of 
M > 0.008 Moyr"^ for the star formation rate averaged over 
the past ~ 2 Myr. From an analysis of the optically visible 
populations of stars in the Carina Nebula, they derived an av- 
erage star formation rate over the past 5 Myr of M 0.010 - 
0.017 Moyr"^ In our analysis of the YSO detected as point- 
like sources in our Herschel data of the entire CNC (i.e., in- 
cluding the Gu m 3 1 region) we estimat ed a star formation rate 
0.017 Moyr"^ (iGaczkowski et al.l l2013l) . which is an average 
over the last few 10^ years. The error associated to the star for - 
mation rate is about 50% (details in IGaczkowski et al.ir2013h . 
Taking the smaller area co vered by the Spitzer stud ies into ac- 
count, it can be shown (see IGaczkowski et"aDl2013l for details) 
that our Herschel-hsised star forma tion rate estimate is in good 
agreement with the rates derived by iPovich et al.l (1201 ll) . 

In Figure \TU\ we add our data for the C NC to the star for - 
mation rate versus cloud mass diagram from iLada et aP (1201 2l) . 
Since the cloud mass as well as star formation rate of the CNC 
are larger than in all other Galactic regions shown in the origi- 
nal plot, the CNC adds an important new data point in the gap 
between the other Galactic regions and the external galaxies. 

Considering the total cloud mass above the Ak > 0.1 mag 
threshold, the data point for the CNC appears to follow the gen- 
eral relation. The cloud mass above the > 0.8 mag threshold, 
however, appears to be considerably lower than the expectation 
according to the general relation. This implies that the star for- 
mation rate per unit mass in dense clouds is higher in the CNC 
compared to the other regions. In other words, the star formation 
process seems to be exceptionally efficient in the CNC. 

Dividing the total mass in dense clouds of 23 000 M© by 
the star formation rate of 0.017 Moyr"^ yields a characteristic 
timescale of 1.3 Myr. If one would simply assume that the to- 
tal mass of the currently present dense clouds will be turned 
into stars at a temporally constant rate, this process would last 
1.3 Myr. Considering the clearly established fact that star for- 
mation in the CNC is going on since (at least) about 5 Myr (as 
witnessed by the large populations of young stars in the opti- 
cally visible populations), this would imply that the star forma- 
tion process in the CNC will soon come to its close and the ob- 
served dense clouds are the last remaining parts of the original 
clouds. 



^ Ak = 0A mag - Ay = 0.9 mag; A^^ = 0.8 mag - Ay = 7.3 mag. 
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However, we believe that it is probably not correct to think 
in terms of a fixed amount of dense clouds that just turns into 
stars over time. Rather, a large fraction of the dense clouds that 
are present today, may have just recently been actively created 
by the feedback of the massive stars. 

Several studies of the youngest stellar populations in the 
CNC came to the conclusion that the currently ongoing star 
formation process is, to a large degree, spatially restricted to 
the edges and surfac es of irradiated c l ouds, e.g. the tips of 
the numerous pillars (iSmith et al.ll2QlQl: lOhlendorf et"aDl2Q12t 
iGaczkowski et al .1120131) . This suggests a causal connection be- 
tween the star formation process and local cloud compression. 
The ongoing feedback continues to compress (moderately dense) 
clouds and thus can create new dense cloud structures that, once 
they get gravitationally unstable, will collapse and form stars. 

In this dynamic picture, where dense cloud structures are 
constantly created, the ratio of dense cloud mass versus the star 
formation rate is naturally lower than in the case of clouds that 
do not experience significant levels of compression from massive 
star feedback. This may explain the position of the CNC in this 
plot. The CNC shows the largest deviation from the SFR versus 
cloud mass relation than all other plotted Galactic star forming 
regions because the level of massive star feedback is particularly 
high. 

It is interesting to note some of the "(ultra-) luminous in- 
frared galaxies" show deviations in the same sense, i.e. SFRs 
that are higher than expected form the dense cloud mass. This is 
consistent with the idea that a high level of feedback is increas- 
ing the efficiency of star formation. 

Since numerous galactic star forming regions have now been 
observed by Herschel, the analysis of these data will soon show 
whether similarly large star forming complexes with comparable 
levels of feedback do also show a similar behavior between cloud 
mass and star formation rate. 

7. Conclusions and Summary 

Our wide-field Herschel SPIRE and R\CS maps allowed us 
to determine the temperatures, surface densities, and the local 
strength of the FUV irradiation for all cloud structures over the 
entire spatial extent 125 pc in diameter) of the Carina Nebula 
complex at a spatial resolution of ^ 0.4 pc for the first time. The 
main results and conclusions of our study can be summarized as 
follows: 

1 . The global temperature structure of the clouds in the CNC 
is clearly dominated by the radiative feedback from the numer- 
ous massive stars. 

2. We present detailed temperature, column density, and 
FUV intensity maps of particularly interesting small-scale cloud 
structures in the CNC. 

3. The comparison of the Herschel FIR maps to a sub-mm 
map of the CNC shows that the cold gas, traced by the sub-mm 
emission, is more spatially concentrated than the warmer gas, 
traced by Herschel. 

4. Considering the masses of individual clumps determined 
from their FIR dust emission, we find that previous mass esti- 
mates that were obtained from CO molecular line observations 
yielded systematically (by factors of ~ 2 - 4) larger values. This 
diff'erence probably results from uncertainties in the absolute cal- 
ibration of the diff'erent mass tracers and from the diff'erent biases 
(and beam sizes) of the observations. 

5. We investigate the total cloud mass budget in the CNC 
at diff'erent column density thresholds and find that the diff'er- 
ent tracers (dust and CO) yield consistent total masses for the 
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Fig. 10. Plot adopted from Fig.2 in lLada etall (120121) . The red points 
represent the position of the star formation rate and the total mass of 
the entire Carina Nebula Complex. The error in the log (M) is 0.1 (i.e. 
-30% of the total mass), and the error in log {SFR) is 0.2 (i.e. -50% 
of the star formation rate). The empty points represent the cloud mass 
computed over an ext inction threshold of A^^ > 0. 1 and the filled point 
over Ak > 0.8 (as in iLada et aLll2012h . three parallel lines correspond 
to /dg = 1-0. 0-1. and 0.01 (defined in Equation [3]) from left to right, 
respectively. 



moderately dense gas (as traced by ^^CO) and the dense gas (as 
traced by ^^CO and C^^O). 

6. The intensity of the FUV irradiation at most cloud sur- 
faces in the CNC is at least about 1000 times stronger than the 
galactic average value (the so-called Habing field). At some lo- 
cations, e.g. at the surface of the massive cloud to the west of 
the stellar cluster Tr 14, the FUV irradiation is more than 10 000 
times the Habing field. 

7. In a region north of the Carina Nebula, we discover a 
"wave"-like pattern in the column-density map. This pattern 
may be related to a flow of hot gas streaming out of the northern 
part of the bipolar superbubble. 

8. We compare the ratio between the total cloud mass and 
the star formation rate in the CNC to the sample of Galactic and 
extragalactic star forming regions originally presented by Lada 
et al. (2012). This shows that the CNC seems to form stars with a 
higher efficiency than other, typically lower-mass. Galactic star 
forming regions. We suggest that this is a consequence of the 
particularly strong radiative feedback in the CNC. The compres- 
sion of the clouds actively and continuously creates dense clouds 
out of lower density diff'use clouds. In this way, the total amount 
of dense clouds is continuously replenished. 

9. We present the first detailed FIR mapping and mass esti- 
mates for two cloud complexes seen at the edges of our Herschel 
maps, which are molecular clouds in the Galactic background at 
distances between 7 and 9 kpc. 
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Table 3. Cloud mass budget (in M©) using different methods. 
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i^CO 
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References. (1) This work; r2) lPreibisch et all (1201 Idl) : G) lPreibisch et all (120121) (Paper I); (^4) lPekruhl etaD(l2013h : (^5) lYonekura et aLl(l2005h . 



Science Center, and the HIFI, PACS and SPIRE consortia. This publication uses 
data acquired with the Atacama Pathfinder Experiment (APEX), which is a col- 
laboration between the Max-Planck-Institut fur Radioastronomie, the European 
Southern Observatory, and the Onsala Space Observatory. V.R. thanks A. Stutz 
for her help with the images convolution, H. Linz and M. Nielbock for their ad- 
vises and useful discussions. We thank C. Lada for the permission to reproduce 
Fig. 10. 



References 

Adams, F. C, HoUenbach, D., Laughlin, G., & Gorti, U. 2004, ApJ, 611, 360 
Aniano, G., Draine, B. T., Gordon, K. D., & Sandstrom, K. 2011, PASP, 123, 
1218 

Bally, J., O'Dell, C. R., & McCaughrean, M. J. 2000, AJ, 119, 2919 

Berne, O., Marcelino, N., & Cemicharo, J. 2010, Nature, 466, 947 

Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 115, 81 

Brooks, K. J., Cox, P, Schneider, N., et al. 2003, A&A, 412, 751 

Burkert, A. & Hartmann, L. 2012 

Caswell, J. L. & Haynes, R. F. 1987, A&A, 171, 261 

Cersosimo, J. C, Mader, S., Figueroa, N. S., et al. 2009, ApJ, 699, 469 

Clarke, C. J. 2007, MNRAS, 376, 1350 

Deharveng, L., Zavagno, A., & Caplan, J. 2005, A&A, 433, 565 
Draine, B. T. & Lee, H. M. 1984, ApJ, 285, 89 
Fang, M., van Boekel, R., King, R. R., et al. 2012, A&A, 539, A119 
Fontani, F., Beltran, M. T., Brand, J., et al. 2005, A&A, 432, 921 
Gaczkowski, B., Preibisch, T., Ratzka, T., et al. 2013, A&A, 549, A67 
Goldsmith, P F., Heyer, M., Narayanan, G., et al. 2008, ApJ, 680, 428 
Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3 
Griffith, M. R. & Wright, A. E. 1993, AJ, 105, 1666 
Gritschneder, M., Burkert, A., Naab, T., & Walch, S. 2010, ApJ, 723, 971 
Gritschneder, M., Naab, T., Walch, S., Burkert, A., & Heitsch, F. 2009, ApJ, 694, 
L26 

Habing, H. J. 1968, Studies of physical conditions in H I regions. 
Hennemann, M., Motte, F., Schneider, N., et al. 2012, A&A, 543, L3 
HoUenbach, D. J., Takahashi, T., & Tielens, A. G. G. M. 1991, ApJ, 377, 192 
Hollenbach, D. J. & Tielens, A. G. G. M. 1999, Reviews of Modem Physics, 71, 
173 

Kramer, C, Cubick, M., Rollig, M., et al. 2008, A&A, 477, 547 

Lada, C. J., Forbrich, J., Lombardi, M., & Alves, J. F. 2012, ApJ, 745, 190 

Minier, V., Tremblin, P, Hill, T., et al. 2013, A&A, 550, A50 

Mottram, J. C, Hoare, M. G., Lumsden, S. L., et al. 2007, A&A, 476, 1019 

O'dell, C. R., Wen, Z., & Hu, X. 1993, ApJ, 410, 696 

Ohlendorf, H., Preibisch, T., Gaczkowski, B., et al. 2012, A&A, 540, A81 

Ohlendorf, H., Preibisch, T., Gaczkowski, B., et al. 2013, A&A, 552, A14 

Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943 

Ott, S. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 434, 
Astronomical Data Analysis Software and Systems XIX, ed. Y. Mizumoto, 
K.-I. Morita, M. Ohishi, 139 

Palmeirim, P, Andre, P, Kirk, J., et al. 2013, A&A, 550, A38 



Pekruhl, S., Preibisch, T., SchuUer, F., & Menten, K. 2013, A&A, 550, A29 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, LI 
Poglitsch, A., Waelkens, C, Geis, N., et al. 2010, A&A, 518, L2 
Povich, M. S., Smith, N., Majewski, S. R., et al. 2011, ApJS, 194, 14 
Preibisch, T., Hodgkin, S., Irwin, M., et al. 2011a, ApJS, 194, 10 
Preibisch, T., Ratzka, T., Gehring, T., et al. 2011b, A&A, 530, A40 
Preibisch, T., Ratzka, T., Kudema, B., et al. 2011c, A&A, 530, A34 
Preibisch, T., Roccatagliata, V., Gaczkowski, B., & Ratzka, T. 2012, A&A, 541, 
A132 (Paper I) 

Preibisch, T., SchuUer, F., Ohlendorf, H., et al. 201 Id, A&A, 525, A92 
Ricci, L., Robberto, M., & Soderblom, D. R. 2008, AJ, 136, 2136 
Richling, S. & Yorke, H. W. 2000, ApJ, 539, 258 

Rivera-Ingraham, A., Martin, P. G., Polychroni, D., et al. 2013, ApJ, in press 

Roccatagliata, V., Bouwman, J., Henning, T., et al. 2011, ApJ, 733, 113 

Roussel, H. 2012, ArXiv 1205.2576 

RusseU, D. 2003, A&A, 397, 133 

RusseU, D. & Castets, A. 2004, A&A, 417, 107 

Schneider, N. & Brooks, K. 2004, PAS A, 21, 290 

Schneider, N., Motte, R, Bontemps, S., et al. 2010, A&A, 518, L83 

Schneider, N., Stutzki, J., Winnewisser, G., Poglitsch, A., & Madden, S. 1998, 

A&A, 338, 262 
Smith, N. 2006, MNRAS, 367, 763 
Smith, N. & Brooks, K. J. 2007, MNRAS, 379, 1279 

Smith, N. & Brooks, K. J. 2008, in Handbook of Star Forming Regions, Volume 

II: The Southern Sky, ed. Reipurth, B., 138 
Smith, N. & Conti, P S. 2008, ApJ, 679, 1467 
Smith, N., Egan, M. P, Carey, S., et al. 2000, ApJ, 532, L145 
Smith, N., Povich, M. S., Whitney, B. A., et al. 2010, MNRAS, 406, 952 
Smith, N., Stassun, K. G., & Bally, J. 2005, AJ, 129, 888 
Stutz, A., Launhardt, R., Linz, H., et al. 2010, A&A, 518, L87 
Townsley, L. K., Broos, P S., Chu, Y.-H., et al. 2011a, ApJS, 194, 15 
Townsley, L. K., Broos, P S., Corcoran, M. R, et al. 2011b, ApJS, 194, 1 
Tremblin, P, Audit, E., Minier, V, Schmidt, W., & Schneider, N. 2012, A&A, 

546, A33 
Vallee, J. P 2008, AJ, 135, 1301 

Yonekura, Y, Asayama, S., Kimura, K., et al. 2005, ApJ, 634, 476 



Appendix A: Cloud complexes at the edges of the 
Herschel maps 

Besides the clouds associated to the Carina Nebula and the 
Gum 31 region, our Herschel maps also show several cloud 
structures at the periphery of the field of view. The two most 
prominent peripheral cloud complexes are seen near the south- 
eastern edge and the western edge of our map. Since no detailed 
FIR observations of these clouds seem to be available so far in 
the literature, we briefly describe the FIR morphology of these 
cloud complexes as revealed by our Herschel maps and discuss 
identifications with observations at other wavelengths. 
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Appendix A.1: The G289.0 - 0.3 cloud complex 

A cloud complex extending over about 35' x 35' is seen around 
the position (Qfj2ooo, ^J20oo) = (10^ 56^ 30\ -60° 06' 00") , near 
the south-eastern edge of our SPIRE maps. Optical images of 
this region show numerous field stars, but no prominent cloud 
features. The cloud complex we see in the Herschel maps is re- 
lat ed to the radio-dete c ted m olecular cloud G289.0-0.3 described 
by Russ eil & CastetsI (12004 ). With a radial velocity of Vrad = 
+21.5 km/s, this cloud complex seems to be clearly unrelated to 
the Carina Nebula cloud complex (which has Vrad ~ -26 km/s). 

Besides the diff'use cloud emission, our Herschel maps re- 
vealed 52 individual p oint-like sources, the po sitions and fluxes 
of which are given in Gacz kowski et al.l ([2013). Several of these 
can be identified with IRAS or MSX point- sources, some of 
which are marked in Fig. lA.ll The brightest compact Her- 
schel source in this area can be identified with the giant HII 
region PMN J1056-6005, whic h was detected in the P arkes- 
MIT-NRAO 4850 MHz survey (lGriffith & Wright) ll 99 3h . and 
se ems to be identical to the H I I region GAL 289.06-00 .36 Hsted 
in lCasweU & HavnesI (Il987h . ICersosimo et al.l (l2009h detected 
emission in the H166Qf radio recombination line from this re- 
gion and determined a distance of Z) = 7.1 ± 0.3 kpc. We note 
that this distance value is consisten t with the pre diction of the 
velocimetric model of the Galaxy by ! Vailed (l2008l) for the galac- 
tic longitude and the measured radial velocity of this molecular 
cloud complex. 

At the southern tip of the cloud complex, the Be star CPD- 
59 2854 IRAS 10538-5958) is seen as a bright compact far- 
infrared source. In the northern part of this complex, one of 
the bright Herschel point-like sources is the extended 2MASS 
source 2MASX J10543287-5939178, which is Hsted as a Galaxy 
in the S IMBAD catalog, but m ay in fact be an ultra-compact H II 
region dBronfman et al.ll 1996b . 

Since this cloud complex is located near the edge of our Her- 
schel maps, and a considerable fraction is only covered by the 
SPIRE maps, but not by PACS, we cannot provide a complete 
characterization of this cloud complex. In our temperature map, 
the dust temperatures in these clouds ranges from about 1 8 K to 
about 25 K. Integrating our column density map (see Sectior[3]) 
over a 16' x 19' area, covering that part of the cloud complex 
that is observed by PACS and SPIRE, we obtain a total mass 
of 65 000 Mq. Due to the incomplete spatial coverage, this is 
clearly a lower limit to the total cloud mass of this complex. 

This rather high cloud mass and the remarkably high number 
of detected point-like Herschel source^ clearly show that this 
cloud complex is a region of active massive star formation that 
is worth to be studied in more detail. 

Appendix A.2: The G286.4- 1.3 cloud complex 

Near the western edge of our PACS maps, a cloud complex 
extending over about 16' x 19' is seen around the position 
(q^j20oo, ^j20oo) = (10^ 34^ 00^ -59° 47' 00"). Optical images 
of this region show several dark clouds. 

The cloud complex seems to be related to the radio-detected 
molec ular cloud G286.4-1.3 discussed by iRusseil & CastetsI 
(|2004|) . The radial velocity of Vrad = +39.9 km/s is clearly 

^ We note that due to the larger distance, the Herschel point- source 
detection limit corresponds to about 10 times larger FIR fluxes than for 
the young stellar objects in the Carina Nebula; this suggests that most 
of the point-like FIR sources detected in this cloud complex must be 
relative massive (^ 10 M©) protostellar objects (or, alternatively, very 
compact and rich clusters of young stellar objects) 



distinct from that of the Carina Nebula cloud complex (Vrad ~ 
-26 km/s). We note that no CO emission from these clouds was 
detected in the survey of Yonekura et al. ( 2005). 

Besides the difl'use cloud emission, our Herschel maps re- 
vealed 15 individual point-like sources in these clouds, the 
positio ns and fluxes of which are given in iGaczkowski et al.l 
(l2013l) . The brightest emission peak in this complex can be 
identified with th e ultra-compact H II region IRAS 10320-5928 
(lBronfmanetal.lll996k . 

The southern part of the cloud complex contains the FIR 
source IRAS 10317-5936, which is identified with the HII re- 
gion GAL 286.42-01.48 and Hsted as part of the star-forming 
complex 286.4-01.4 in Russeil (2003). It contains the two can- 
didate massive young s tellar objects 2 86.3938-01.3514 1 and 2 
(lMottrametal.li2007h . iFontani et aP (i2005) determined a dis- 
tance of Z) = 8.88 kpc to this cloud. We note that this distance 
value is consiste nt with the prediction of the velocimetric model 
of the Galaxy by I Vailed (12008 ) for the galactic longitude and the 
measured radial velocity of this molecular cloud complex. 

In the north-western part of this complex, another accumula- 
tion of compact clouds is found, two of which can be identified 
with IRAS sources. 

The main part, but not the full extent of this cloud complex is 
covered by our PACS and SPIRE maps. The dust temperatures 
range from about 20 K to ^ 35 K. Integrating our column density 
map yields a total mass of ~ 93 000 M© within a 16' x 19' area, 
if we assume the above mentioned distance of 8.88 kpc. This is 
a strict lower limit, since our data cover only some part of the 
cloud complex. 

This rather high cloud mass and the presence of numerous 
bright point-like Herschel sources (with fluxes up to 530 Jy, 
corresponding to FIR luminosities of up to ~ 25 000 L©) shows 
that this complex is a site of active massive star formation and 
certainly worth to be studied in more detail. 
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